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Abstract. Let m,n > 2 be positive integers. Denote by M m the set ofmxm complex matrices 
and by w(X) the numerical radius of a square matrix X. Motivated by the study of operations on 
^ . bipartite systems of quantum states, we show that a linear map </> : M mn — ► Mmn satisfies 

Q . w(<p(A ® B)) = w(A ® B) for all A e M m and B £ M n 

C^) ' if and only if there is a unitary matrix U £ M mn and a complex unit £ such that 

(t>(A ®B)= £U((pi{A) ® <p 2 (B))U* for all ^ G M m and B <E M n , 

where tpk is the identity map or the transposition map X M- X for k = 1,2, and the maps 
rj^ and <y32 will be of the same type if m, n > 3. In particular, if m, n > 3, the map corresponds to 

an evolution of a closed quantum system (under a fixed unitary operator), possibly followed by a 
transposition. The results are extended to multipartite systems. 
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1. Introduction and preliminaries 

\ Let M n be the set of n x n complex matrices for any positive integer n. For A € M n , define (and 

| denote) its numerical range and numerical radius by 



W(A) = I u*Au : u G C n , u*u = 1 } and w(A) = sup{|^| : ^ £ W(A)}, 



respectively. The study of numerical range and numerical radius have a long history and is still 
under active research. Moreover, there are many generalizations motivated by pure and applied 
topics; see [UEJE]. 

By the convexity of the numerical range, 

W{A) = {tr (Auu*) : u G C n , u*u = 1} = {tr (AX) : X G D n }, 

where D n is the set of density matrices (positive semidefinite matrices with trace one) in M n . In 
particular, in the study of quantum physics, if A G M n is Hermitian corresponding to an observable 
and if quantum states are represented as density matrices, then H^(^4) is the set of all possible 
measurements under the observables and w(A) is a bound for the measurement. If A = Ax + 1A2, 
where A\,Ai G M n are Hermitian, then VF(^4) is the set of the joint measurement of quantum 
states under the two observables corresponding to A\ and A%. 
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Suppose m, n > 2 are positive integers. Denote by A ® B the tensor (Kronecker) product of the 
matrices A G M m and -B G M n . If 4 and B are observables of two quantum systems, then A®B is 
an observable of the composite bipartite system. Of course, a general observable on the composite 
system corresponds to C G M mn , and observable of the form A® B with 4 G M m , B G M n is a 
very small (measure zero) set. Nevertheless, one may be able to extract useful information about 
the bipartite system by focusing on the set of tensor product matrices. In particular, in the study 
of linear operators <\> : M mn — > M mn on bipartite systems, the structure of 4> can be determined by 
studying cf>(A <g> B) with A G M m , B G M n ; see [H El 02 E| and their references. 

In this paper, we determine the structure of linear maps <j> : M mn — > M mn satisfying w(A<g>B) = 
w{(j){A £5 B)) for all A G M m and B G M n . We show that for such a map there is a unitary matrix 
U G M mn and a complex unit £ such that 

0(4 ® J3) = £E%i(A) (g) if 2 {B))U* for all A G M m and 5 G M n , 

where is the identity map or the transposition map X \— > X 1 for k = 1,2, and the maps (pi 
and (f2 will be of the same type if m, n > 3. In particular, if m, n > 3, the map corresponds to 
an evolution of a closed quantum system (under a fixed unitary operator), possibly followed by a 
transposition. 

The study of linear maps on matrices or operators with some special properties are known as 
preserver problems; for example, see [9] and its references. In connection to preserver problems on 
bipartite quantum systems, it is quite common that if one considers a linear map eft : M mn — > M mn 
and imposes conditions on <f>{A ® B) for A G M m , B G M n , then the partial transpose maps 
A ® B i— > A IS) B and A ® B \— > A 1 ® B are admissible preservers. It is interesting to note that 
for numerical radius preservers and numerical range preservers in our study, if m, n > 3, then the 
partial transpose maps are not allowed and that the (linear) numerical radius preserver on M mn 
will be of the standard form 

X ^ £V*XV or X ^ ^riV 

for some complex unit £ and unitary V G M mn . This is the first example of such results in this line 
of study. It would be interesting to explore more matrix invariant or quantum properties that the 
structure of preservers on Af mn can be completely determined by the behavior of the map on the 
small class of matrices of the form A ® B G M mn with A G M m , B G M n . 
In our study, we also determine the linear map 4> : M mn — > M mn such that 

W{A ®B) = W((p(A ® B)) for all (A, B) G M m x M n . 

We will denote by X t the transpose of a matrix X G M n and X* the conjugate transpose of a 

(n) 

matrix X G M n . The n x n identity matrix will be denoted by I n . Let E\- G M n be the matrix 

whose (z,j)-entry is equal to one and all the others are equal to zero. We simply write Eij = 
if the size of the matrix is clear. 

We will prove our main result on bipartite systems in Section 2 and extend the results to multi- 
partite systems in Section 3. 
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2. Bipartite systems 
The following example is useful in our discussion. 

/0 2 0\ 

Example 2.1. Suppose m, n > 3. Let A = X © O m _ 3 and £ = X © O n _ 3 with X = 1 J . 

\0 0/ 

Then ^4 © £? is unitarily similar to 

(o i) e (jj ° v 2 )©o mn _ 7 , 

and ^4 (g) £>* is unitarily similar to 

One readily checks (see also [10]) that <g> B) and <g> £?*) = W(A t <g) B) are circular disks 
centered at the origin with radii w(A<S> B) and w(A<g> B l ), respectively. Moreover, we have 

w{A®B) = A max (^ ® B + (A<8> Bf) /2 = v^o" 

> 2.0000 = \ max (A®B t + (A ®B t ) t )/2 = w{A®B t ). 



In the following, we first determine the structure of linear preservers of numerical range using 
the above example and the results in [2]. 

Theorem 2.2. The following are equivalent for a linear map <p ■ M mn — > M mn . 

(a) W((j){A ® B)) = W{A ® B) for any A G M m and B G M n . 

(b) There is a unitary matrix U G M mn such that 

<p(A ®B) = U(<pi(A) <g> ip 2 (B))U* for all A G M m and B G M n , 

where <pk is the identity map or the transposition map X i— > X 1 for k = 1,2, and the maps 
(pi and (f2 will be of the same type if m,n > 3. 

Proof. Suppose (b) holds. If m,n > 3, then the map has the form C h-> UCU* or C H> UC 1 !!*. 
Thus, the condition (a) holds. If m = 2, then ^4* and A are unitarily similar for every A G M2. So, 
W{A ®B) = W(A l ® 5) for any 5 G M n . Hence, the condition (a) holds. Similarly, if n = 2, then 
(a) holds. 

Conversely, suppose that W(0(^4 ® 5)) = W(A ® B) for all A G M m and B G M n . Assume for 
the moment that A ® B G M mn is a Hermitian matrix. Then 

This yields that 0(^4 (g> -B) is a Hermitian matrix, as well. Thus, 4> maps Hermitian matrices 
to Hermitian matrices and preserves numerical radius, which is equivalent to spectral radius for 
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Hermitian matrices. By Theorem 3.3 in [2], we conclude that eft has the asserted form on Hermitian 
matrices and, hence, on all matrices in M mn . However, if m,n > 3, then, by Example 12.11 neither 
the map A <g> B \— > A <g> B l nor the map A ® B i— > A 1 ® B will preserve the numerical range. So, the 
last statement about (pi and </?2 holds. □ 

Next, we turn to linear preservers of the numerical radius. We need the following (well-known) 
lemma to prove our result. We include a short proof for the sake of the completeness. 

Lemma 2.3. Let A G M n with w(A) = |x*^4x| = 1 for some unit x G C n . Then for any unitary 
U G M n with x being its first column, there exists some y G C n_1 such that 



(1) U*AU = x*Ax 



1 y* 
-y * 



Proof. Write (x* Ax)" 1 A = G + iH with G and H Hermitian. Then the largest eigenvalue of G is 
1 with x as its corresponding eigenvector and U*GU = [1] © G\ for some Gi G M n -\. Moreover, 
the (l,l)-entry of iU*HU is since w(A) = 1. Since U*HU is a Hermitian matrix we have 

iU*HU = | ° y ). Thus, U*AU has the claimed form. □ 

-y * 

Theorem 2.4. The following are equivalent for a linear map (j) : M mn — > M mn . 

(a) w((p(A (g) B)) = w; (j4 ® B) /or any A £ M m and B G M n . 

(b) There is a unitary matrix U G M mn and a complex unit £ such that 

4>(A ®B)= £U((pi(A) ® ip 2 (B))U* for all A G M m and 5 G Af n , 

where ifk is the identity map or the transposition map X i— >■ X /or A; = 1,2, and t/ie maps 
tjCi and cp2 will be of the same type if m,n > 3. 

Proof. The implication (b) =>■ (a) can be verified readily. Now, suppose (a) holds and let Bij = 
4>{Ea (g) Ejj) for 1 < i < m, 1 < j < n. According to the assumptions, for all 1 < % < m, 1 < j < n, 
there is a unit vector Uij G C mn and a complex unit £ij such that u*jBijUij = We will first 
show that there exists a unitary matrix U G M mn and complex units £y such that 

B^ = £ijU(Eu (g) Ejj)U*, 1 < i < m,l < j < n. 

We divide the proof into several claims. 

Claim 1. Suppose (i, j) 7^ (r, s). If n G C mn is a unit vector such that \u*Biju\ = 1, then i? rs n = 0. 



. _ J8 T of C — „,* R 

*J 1 



Proof. Suppose u is a unit vector such that u*BijU = e ld . Let £ = u*B rs u. We first consider the 
case if z = r. For any [i G C with |//| < 1, 

(2) 1 = w(E u (g) (Ejj + ui? ss )) = + ufi rs ) > \u*{Bij + nB rs )u\ = \e ie + ^\. 

Then we must have £ = 0. Otherwise, \e lS + /Lt£| = |e i61 + e i6, |£|| > 1 if one chooses u = e* e £/|£|. 
Suppose £ = 0. Then the inequality in (|2|) become equality and w(Bij+[iB rs ) = \u*(Bij+fiB rs )u\ = 
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1. By Lemma |2, 3\ there is € C mn -! such that 

U*B i:j U + fiU*B rs U = U*(B lj + nB rs )U = e ie ( 1 y " ) , 

\-Vfi * J 

where U is a unitary matrix with u as its first column. As the above equation holds for any /16 C 
with \fi\ < 1, U*B rs U must have the form ^jj ^ J , or equivalently, i? rs u = 0. Similarly, we can 
prove the case if j = s. 

Now we consider the case when i ^ r and j ^ s. By the previous argument, we have Bi s u = 
B r jU = 0. Then for any \i 6 C with < 1, 

1 = w{{Eu + E rr ) ® (Ejj + nE„)) 
= w(Bij + B rj + n(B is + B rs )) 
> \u*(B i:j + Brj + n{B rs + B rs ))u\ 
= \u* Biju + nu* B rs u\ 
= \e i6 + rt\. 

It follows that ^ = u*B rs u = and hence w(Bij + 5 r j + fi(Bi s + B rs )) = \u*Biju\ = 1. By Lemma 
12.31 we conclude that (Sj S + B rs )u = and thus, B rs u = 0. □ 



Claim 2. Suppose / (r, s). If Uij,u rs £ C mn are two unit vectors such that \u*jBijUij\ = 
\u* s B rs u rs \ = 1, then u*-u rs = 0. 

Proof. Suppose u rs = auij+/3v with a = u* s Uij and f3 = u* s v, where v is a unit vector orthogonal 
to Uij. Notice that \a\ 2 + |/3| 2 = 1. By Claim 1, u*jB rs = and B rs uij = and so 

1 = \u* rs B rs u rs \ = \P\ 2 \v*B rs v\ < \p\ 2 w(B rs ) = |/3| 2 < 1. 
Thus, |/3 1 = 1 and hence u* s Uij = a = 0. □ 

Claim 3. Let U = [u\\ ■■■ u\ n Ui\ • • ■ U2 n •■• u m i ••• Umn]- Then U*U = I mn and U*BijU = 
£,ij(Eu ® Ejj) for all 1 < i < m and 1 < j < n. 

Proof. By Claim 2, {u^ : 1 < i < m, 1 < j < n} forms an orthonormal basis. Thus, U*U = I mn . 
Next by Claim 1, u* s BijUke = for all (r,s) and (k,£), except the case when (r, s) = (k,£) = (i,j). 
Therefore, the result follows. □ 

According to our assumptions and by Claims 1, 2, 3, we see that up to some unitary similarity 

<f>(Eu ® Ejj) = iij{Eu ® Ejj) 

for 1 < i < m,l < j < n and some complex units Now, for any unitary X € M m , using the 
same arguments as above, there exists some unitary Ux and some complex units fiij such that 

<f>(XE u X* ® Ejj) = iHjUxiEa ® E 3j )U x 
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for all 1 < i < m and 1 < j < n. We see that 4>{XEnX* ®Ejj) is a rank one matrix with numerical 
radius one. If 7 > 0, then 

w(<f>((XEuX* + jl m ) ® Ejj)) = 1 + 7- 

Thus, <p(XEaX* ® Ejj) has the form Ri t x ® Ejj for some Ri.x G M m . Since this is true for any 
1 < i < m and unitary X G M m , we have 

0(^4 <g> Ejj) = (fj(A) (£> 

for all matrices A G M m and some linear map ipj. Clearly, ipj preserves numerical radius and, 
hence, has the form 

A ^ ij WjAW* or A h> Wj A 1 W* 

for some complex unit £j and unitary Wj G M m . In particular, ipj{I m ) = £j/ m and (p(I mn ) = I m ®D 
for some diagonal matrix D G M n . Using the same arguments as above, we can show that 

4>{E ii ®B) = Eu®<Pi(B) 

for all matrices B G M n and some linear map ipi of the form 

5 ^^WiBW* or S ^ eiW-fl*^*, 

where £j is a complex unit and Wi G M n a unitary matrix. Therefore, we have (fi(I n ) = (,il n and 
4>(Imn) = D®I n fo r some diagonal matrix D G M m . Since 4>(I m n) = D®I n = I m ®D, we conclude 
that 4>(I mn ) = C-^mn for some complex unit £. For the sake of the simplicity, let us assume that 
<t>{Im.n) = Imn- Then (j){Eu <g> Ejj) = En <g> Ejj for all 1 < i < m, 1 < j < n. 

For any Hermitian matrices A G Af m and B G M n , suppose their spectral decompositions 
are A = XD\X* and B = YD2Y*. Repeating the above argument and using the assumption 
4>(Imn) = Imn, one sees that there exists a unitary matrix Ux,y such that 

(j)(XEuX ® YEjjY*) = Ux,y(XEuX* <g> KE^Y*)^^, 1 < i < m, 1 < j < n, 

and, hence, (/>(t4 ® B) = Ux,y{A <g> B)U XY . So, maps Hermitian matrices to Hermitian matrices 
and preserves numerical range on the tensor product of Hermitian matrices. Thus, by the same 
argument as in the proof of Theorem 12.21 <p has the asserted form on Hermitian matrices and, 
hence, on all matrices in M mn . If m,n > 3, we can use Example 12.11 to conclude that tp\ and 922 
should both be the identity map, or both be the transpose map. The proof is completed. □ 

3. Multipartite systems 

In this section we extend Theorem 1 2 . 2 1 and Theorem l2.4l to multipartite systems M ni &>- ■ •®M nm , 
m > 2. 

Theorem 3.1. Let m,...,n m > 2 be positive integers and N = Jlj=i n j- The following are 
equivalent for a linear map <f> : Mjy — > Mn ■ 

(a) <g> ■ ■ ■ <g> A m )) = W{At ® ■ ■ ■ ® A m ) for any (A u . . . , A m ) G M ni x • • • x M nm . 
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(b) There is a unitary matrix U G such that 

(3) <£(Ai ® ■ • ■ ® 4*) = U(tp 1 (A 1 )®---®<p rn (A m ))U* 

for all (A\, . . . , A m ) G M ni x • • • x M„ m , where (p^ is the identity map or the transposition 
map X i— > X 1 for k = 1, . . . , m, and the maps cpj are of the same type for those j 's such 
that rij > 3. 

Proof. The sufficient part is clear. For the converse, as in the proof of Theorem 12.21 consider 
Hermitian matrix A = A\ ® ■ ■ ■ ® A m with Aj G H n . for j = 1, . . . , m. By [2j Theorem 3.4], <p has 
the asserted form on Hermitian matrices and, hence, on all matrices in M ni <S> • • • M nm . 

However, if ni,nj > 3 with i < j, let Ai = X © O ni s and Aj = X © O nj _3, where X is defined 
as in Example 12. 1\ and A = £n G M nfc for A; 7^ i,j. Then w(A\ <8> • • • <8> -Am) = V4.25 and 
(8) ■ ■ ■ ® ® A* ® A+l ® ' ' ' ® -4m) = 2.0000. Thus, 

® • • • ® A m ) / W(Ai ® • • • <g> Aj-x ® A*- A+i ® ■ ■ ■ ® 
and we see that the last statement about (p^ holds. □ 

Theorem 3.2. Let ni,...,n m y 2 be positive integers and N — 11^=1 The following are 
equivalent for a linear map (f> : Mjv — > M^. 

(a) w(cp(A 1 • • • ® A m )) = w{Ai <g> ■ ■ ■ ® A m ) for any (A\, . . . , A m ) G M ni x • • • x M nm . 

(b) There is a unitary matrix U G Mn and a complex unit £ smc/i £/ia£ 

(4) ® ■ ■ ■ ® A m ) = Z,U(ipi(M) <g> • • • <g> <p m (A m ))U* 

for all (A\, . . . ,A m ) G M ni x • • • x M nm , where ip^ is the identity map or the transposition 
map X I—)- X t for k = 1, . . . , m, and the maps cpj are of the same type for those j 's such 
that nj > 3. 

Proof. The sufficiency part is clear. We verify the necessity. First of all, one can use similar 
arguments as in the proof of Theorem 12.41 (see Claim 1, Claim 2, and Claim 3) to show that up to 
some unitary similarity, 

<8> E i2i2 tg> ••• ® E imim ) = Uii 2 -im( E iih ® E ii%i ® ••• ® #w m ) 
for all 1 < if. < with 1 < k < m and some complex units ii 1 i 2 ---i m - Below we give the details of 
the proof for the case M ni <g> M n2 ® M„ 3 . One readily extends the arguments to the general case. 

Denote Bijk = (f>(Ea <S> Ejj ® E^k) for 1 < i < m, 1 < j < ri2, 1 < < 713 and let iV = n\U2n^. 
According to the assumptions, for any 1 < i < n\,l < j < n2,l < k < n^, there is a unit vector 
Ujjfc G and a complex unit such that u*- k Bij^Uijk = 

Claim 4. Suppose (i,j,k) ^ (r,s,t). If n G is a unit vector such that \u*Bij^u\ = 1, then 
B rst u = 0. 

Proof. Suppose u is a unit vector such that u*Bijf~u = e %6 '. Let £ = u*B rs tu. First, assume that 
<5j r + 5j S + 5 kt = 2, where 5 a b equals to 1 when a = b and zero otherwise. In other words, we assume 
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that exactly two of the three sets {i,r},{j,s},{k,t} are singletons. Without loss of generality, 
assume that i = r, j = s, and k 7^ t. For any fj, 6 C with |/x| < 1, 

1 = w(B ijk + nB rst ) > \u*(B ijk + nB rst )u\ = \e %e + fj£\. 

Then we must have £ = 0. Furthermore, wlB^ + ^.B rst ) = \u*{Bij k + /j,B rst )u\ = 1. By Lemma 

12.31 one conclude that U*B rs tU has the form ( n ], where U is a unitary matrix with u as its 

v° *J 

first column, and hence B rs tu = 0. 

Next, suppose 5i r + 5j s + S^t = 1, say i = r. By the previous case, B isk u = B^u = 0. Then for 
any fi £ C with < 1, 

1 = w(Bijk + B is k + ^{Bijt + Bist)) 
> |u*(.B ijfc + 5 isfc + /i(fiij t + B ist ))u\ 
= \u*B ijk u + p,u*B ist u\ = \e td + 

It follows that f; = u*Bi s tu = and hence w(Bijk + I?j s fc + fi(Bijt + Bi s t)) = \u*Bijku\ = 1. By 
Lemma |2T3| we conclude that (Bij t + B ist )u = and thus, -Bj^u = 0. 

Finally, suppose <5j r + 5j s + <5fct = 0. By the previous cases, B^u = B r j k u = B rs kU = B^u = 
Bistu = B r jtu = 0. Then for any [i G C with < 1, 

1 = w(Bijk + B is k + B r jk + B rs k + n{Bijt + Bist + B r j t + B rs t)) 
> \u*(Bij k + Bisk + Brjk + B rs k + n{Bijt + Bist + B r j t + B rs t))u\ 
= \u*Bijt~u + fiu*B rs tu\ = \e ie + 

Then by a similar argument, one conclude that £ = and B rs tu = 0. □ 
Claim 5. Suppose (i,j,k) 7= (r,s,t). If Uijk,u rs t G C mn are two unit vectors such that 

\ U ijk B ijk u ijk\ = \u*st B rstU rs t\ = 1, then U* jk U rs t = 0. 

Proof. Suppose u TS t = cxuijk + (3v with and (3 = u* st v, where v is a unit vector 

orthogonal to Uijk- Note that \a\ 2 + |/3| 2 = 1. By Claim 1, u*- k B rs t = and B rs tUijk = and so 

1 = \u* rst B rst u rst \ = \/3\ 2 \v*B rst v\ < \f3\ 2 w(B rst ) = \/3\ 2 < 1. 
Thus, |/3| = 1 and hence u* st Uijk = a = 0. □ 

Claim 6. Let 

U=[llm ••• Un n3 Ui21 ••• Ul2n 3 ■"• ^ln 2 n 3 '•• ^211 -- - U2n 2 n 3 U ni n ••• U nin2ns \. 

Then U*U = In and U*Bi jk U = Cijk(Eu ® E'jj ® -Efefc) for all 1 < i < m, 1 < j < ri2, 1 < k < n 3 . 

Proof. By Claim 2, {u^ : 1 < i < n\, 1 < j < 112, 1 < k < 113} forms an orthonormal basis 
and thus U*U = In- Next by Claim 1, u* st Byj.^ p? = for all (r,s,t) and (£,p,q), except the case 
when (r,s,t) = (£,p,q) = (i,j,k). Therefore, the result follows. □ 
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Similarly, for any unitary X G M ni , there exists some unitary Ux and some complex units 
fJ>iii 2 ~i m such that 

ftXEi^X* ® E i2i2 ® • • • ® E imi J = fi lli2 ... lm U x (XE hll X* ® £^ 2 ® • • • ® E lmlm )U x 

for all 1 < ifc < nfc with 1 < fe < m. We see that (^(XEi^X* ® -Ej 2 j 2 ® ■ ■ ■ ® Ei m i m ) is a rank one 
matrix with numerical radius one. If 7 > 0, then 

^((JO^X* + 7l ni ) ® ^ 2i2 ® • • • ® £ im iJ) = 1 + 7. 

Thus, 4>{XE ixil X* ® Ej 2 j 2 ® • • • ® E imim ) has the form R®E i2i2 ® • • • ® E imim ) for some i? G M ni . 
Since this is true for any unitary X £ M ni , we have 

0(4 ® E i2i2 ® • • • ® E im i m ) = <p i2 ... im (A) ® #i 2 i 2 ® • • • ® £; mim 

for all Hermitian matrices A G M ni and some linear map (fi 2 ---i m - Clearly, <fi 2 ---i m preserves numer- 
ical radius and, hence, has the form 

,1 • ► W i2 ...i AW? , or ,1 • ► U , ,1'U V, , 

for some complex unit & 2 ...i m and unitary W i2 ... im G M ni . In particular, (p i2 ... im (I ni ) = & 2 ...i m I m 
and 4>(In) = In x ® D\ for some diagonal matrix D\ G M n2 ...„ m . 

Given 2 < k < m and using the same arguments as above, one can show that 4>(In) = Dkl ® 
In k ® D k2 for some diagonal matrix D kl G M„ 1 ... nfe _ 1 and D fc2 G A/„ fc+1 ...„ m . Since 

4>(In) = Im ®Di = D kl ® J n)fe ® L> fc2 for k = 2, . . . , m, 

we conclude that c/>(In) — f° r some complex unit £. For the sake of the simplicity, let us 
assume that 4>(In) = In- Then 

4>{E ixil ® £ J2 ; 2 ® • • • ® S imim ) = E hil ® E i2i2 ® • • • ® £ imim 

for all 1 < ifc < n k with 1 < k < m. 

For any Hermitian matrix A\ ® ■ ■ ■ ® 4 m G M^r, suppose its spectral decomposition is 

At ® • • • ® A m = XtDtXl ® • • • ® X m L» m X^. 

Repeating the above argument and using the assumption 4>(In) = In, we can conclude that there 
exists a unitary matrix Uxi_ ...x m such that 

(t>{X x Ei lix XZ ® • • • ® X m E imim X* m ) = U Xl ,...,x m {XiE ilix Xl ® • • • ® I m £ !m!m i;)^ lr .. >jfra 

for all 1 < ifc < nfc with 1 < k < m. By linearity, we have 

0(4! ® • • • ® A m ) = U Xl) ...,X m ( A l ® • • • ® 4 m )^ lv .. jXm . 

So, maps Hermitian matrices to Hermitian matrices and preserves numerical range on the tensor 
product of Hermitian matrices. By Theorem 13.11 it has the asserted form on Hermitian matrices 
and, hence, on all matrices in M ni ... nm . If ni,rij > 3, we observe the matrices A{ = X © n ._3 
and Aj = X © nj _3, where X is defined as in Example 12. 1( and A k = En G M nk for k ^ i,j, 
to conclude that (pi and ifj should both be the identity map, or both be the transpose map. The 
proof is completed. □ 
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